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Introduction 


Sound  Speed 
(km/s) 


Figure  1:  Arctic  Ocean  sound-speed  profile  and  ray  diagram. 


The  Arctic  Ocean  is  a  unique  acoustic  environment,  principally  because 
the  sound  speed  increases  monotonicelly  with  depth  and  the  surface  is 
generally  covered  by  an  ice  layer.  Upward  refraction  and  scattering  at  the 
rough  ice-water  interface  produce  reverberation  and  attenuation,  both  of 
which  place  severe  limitations  on  SONAR  performance. 

Figure  1  shows  a  typical  sound-speed  profile  end  ray  diagram.  In  deeper 
waters,  the  maximum  RSR  cycle  distance  can  be  of  the  order  of  50  km.  In 
this  range,  reverberation  and  absorption  in  the  medium  con  be  the  most 
important  factors.  At  longer  ranges,  most  of  the  energy  is  concentrated  in 
rays  with  small  grazing  angles  and  multiple  reflections  at  the  ice-water 
interface  ore  then  involved.  Scattering  causes  losses  in  both  energy  and 
signal  coherence  and  these  losses  con  exceed  absorption  loss  by  more  than 
on  order  of  magnitude. 

The  main  purpose  of  this  report  is  to  e- tend  the  current  global-model  for 
sound  absorption  in  sea  water  (1,21  to  induce  Arctic  regions.  The  model  is 
based  on  the  three  principal  chemical  relaxations  contributing  to  sound 
absorption.  An  approximate  formula  has  been  developed  and  predictions 
compared  favorably  with  experiments  data  from  many  parts  of  the  World 
Ocean.  The  environmental  factors  involved  in  absorption  prediction  ore  pH 
and  temperature.  Investigation  of  the  effects  of  pH  on  propagation  loss 
will  be  a  principal  concern.  However,  reverberation  and  attenuation  due  to 
scattering  will  also  be  addressed  because  of  their  great  importance  in  the 
overall  problem 


Background 


By  the  end  of  WWII,  it  was  known  that  absorption  at  SONAR  frequencies 
is  more  than  order  of  magnitude  greater  in  sea  water  than  in  fresh  water. 
Subsequent  propagation  experiments  provided  accurate  estimates  of  the 
magnitude.  By  the  early  1950’s,  laboratory  resonator  experiments  had 
identified  the  cause  os  a  magnesium  sulfate  relaxation  and  details  of  the 
mechanism  and  the  relaxation  parameters  had  been  worked  out.  In  1962, 
Schulkin  and  Marsh  [31  proposed  a  formula  for  see-water  absorption,  based 
on  both  the  measured  relaxation  parameters  and  field  experiments  in  the 
frequency  range  2-30  kHz. 

In  1965,  Thorp  [4]  reported  sound-channel  propagation  experiments  in 
the  Bermudn-Eleuthere  area  indicating  an  anomaly  at  lower  frequencies. 
The  extra  loss  was  fitted  by  adding  a  I  kHz  relaxation  to  the  S&M  formula 
and  the  result  became  known  as  the  Thorp  formula-  [51. 

Mediterranean  experiments  were  reported  by  Leroy  [61,  showing  a  similar 
anomaly  but  with  somewhat  higher  magnitude  and  relaxation  frequency. 
Experiments  carried  over  the  next  two  decades  in  other  areas  confirmed 
the  high  degree  of  variability  of  the  extra  loss  and  regional  dependence  has 
therefore  become  much  more  critical  [7).  The  major  factor  involved  in  this 
variability  has  been  Identified  as  the  pH  value  [8], 

The  principal  chemical  relaxation  responsible  for  the  pH-dependent  loss 
has  been  shown  to  involve  boric  acid  191.  The  chemical  mechanism  has  been 
identified  ns  the  boric  acid/cerbonate  equilibrium  and  the  parameters  have 
been  measured  in  the  laboratory  by  means  of  the  resonator  method  [10). 

The  laboratory  investigations  also  revealed  a  pH-dependent  relaxation 
involving  magnesium  carbonate,  which  has  been  found  to  ploy  a  minor  but 
significant  role  in  sea  water  absorption  111). 

An  absorption  formula,  based  solely  on  known  chemical  processes,  would 
be  for  too  complex  and  the  accuracy  would  be  limited  os  well.  However, 
since  the  range  of  environmental  parameters  in  the  World  Ocean  is  very 
limited,  simplifying  approximations  can  be  made;  namely,  that  the  losses 
for  both  pH-dependent  relaxations  increase  exponentially  with  pH  and  the 
relaxation  frequencies  increase  exponentially  with  temperature. Thorps 
formula  can  then  simply  be  modified  by  adding  the  third  relaxation  ond 
including  the  required  pH  end  temperature  corrections. 

A  3-relaxotion  formula  has  already  been  developed.  Predictions  based  on 
archival  pH  data  have  been  tested  against  ell  the  available  sound-channel 
with  good  results  1 1,121 


Absorption  Model 


A=A,(MgS04)+A2(B(OH)3)+A3(MgC03) 
An=(s/35)anf2fn/(f2+f„  ) 


a^o.Sxio"0^1"^20 

(pH- 8) 

a2=o.ixio 

,  T/60 

fj  =  50x10 

f  2  =  0.9x  f  0T/7° 

a3=  0.03xlO(pH’8) 

f3  =  4.5x10T/3° 

Atlantic  4*c  pH  8.0 

A=0.007f  2  *0.1  f2/(l*f2)+0.18f2/(62+f2) 

N.Pacific  4*C  pH  7.7 

A=0.007f2*0  05f2/(1+  f2)*0  09  f2/(62+f2) 
Mediterranean  14*C  pH  8.3 
A=0.006f2+0.26f2/(1.42*f2)+0.78f2/(122*f2) 

Red  Sea  22*C  pH  8.2 

A~0.004f2*0.27f2/(182»y2)*1.1f2/(242*f2) 

sub-Arctic  -1*C  pH  8.3 

A=0.01f2*0.17  f2/(C.852*f2)*0.24  f2/(42*f2  ) 

Figure  2:  Simplified  absorption  formulae. 

In  the  3-relaxetion  formula  of  Figure  2,  A  is  in  dB/km,  frequency  f  and 
relaxation  frequencies  fn  are  in  kHz,  temperature  T  is  in  °C  and  pH=B.O  the 
reference  value.  PH  values  in  the  World  Ocean  very  roughly  from  7.7  to  6  3, 
which  corresponds  to  an  absorption  ratio  of  as  much  as  4/1  at  the  lower 
frequencies. 

The  magnesium  sulfate  term  includes  the  depth  factor  D(km),  which  is 
adapted  from  th8  pressure  correction  of  Fisher  and  Simmons  {131.  Depth 
dependencies  of  the  other  two  relaxations  are  not  yet  known,  however, 
measurements  in  both  deep  and  shallow  channels  indicate  that  bone  acid 
effects  are  negligible.  Magnesium  carbonate  effects  may  be  greeter  but 
can  be  neglected  because  its  contribution  is  so  small  Salinity  dependence 
has  been  approximated  os  S/35  where  S  is  parts  per  thousand. 

Specific  coefficients  for  several  experimental  areas,  are  shown  in  the 
bottom  box.  Note  that  the  magnesium  sulfate  terms  are  approximations 
valid  only  for  frequencies  less  than  10  kH2. 


Model  end  date  comparison 


Figure  3:  Thorp's  data  8nd  3-relaxation  model. 

Figure  3  compares  Thorp's  data  with  the  3-relaxetion  model.  Individual 
components  ore  identified  end  the  top  curve  is  their  sum.  The  overall  fit  to 
the  data  with  the  new  model  is  as  good  or  better  than  that  with  the  Thorp 
formula:  A=  40f2/(410(M2)  ♦  0  1f2/(W2)  dQ/kyd. 

Note  that  the  A2  (boric  acid)  coefficient  Is  10$  lower  than  that  of  the 
second  term  in  the  Thorp  formula;  i.e.  the  value  in  dB/km  becomes  equal  to 
the  Thorp  value  In  dB/kyd.  The  differences  in  total  absorption  at  the  lower 
frequencies  are  made  up  by  the  A3  (magnesium  carbonate)  component. 

The  parameter  adjustment  is  mainly  justified  on  the  basis  of  dota-fit, 
the  third  component  being  essential  to  the  model.  When  the  sea-water 
resonator  data  were  fitted  with  a  2-relaxation  model,  there  were  serious 
discrepancies  at  higher  pH  values.  For  example,  at  pH=8.5,  formula  values 
that  were  too  low  by  a  factor  of  more  th8n  two  at  the  lower  frequencies, 
which  was  clear  evidence  of  the  existence  of  a  third  component.  Sea  water 
synthesis  experiments  were  carried  out  and  the  mechanism  was  identified 
as  the  magnesium-carbonate  relaxation.  The  relaxation  parameters  and  the 
temperature  and  pH  dependencies  were  determined  by  measurement. 
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Figure  4:  Model  and  data  comparison. 

Figure  4  compares  the  3-relaaatlon  model  predictions  with  data  from 
sound-channel  experiments  In  the  North  Atlantic  (Thorp),  North  Pacific, 
Mediterranean  Sea  and  estimates  for  the  Arctic  (dashed  curve). 

In  the  North  Pacific  case  [14],  the  lower  value  pHz7.7  reduces  both  the 
Doric  acid  (A2)  and  the  magnesium  cerbonate  (A3)  coefficients  by  a  factor 
of  two  compared  to  the  N.  Atlantic.  Relaxation  frequencies  depend  only  on 
temperature  and  remain  the  same. 

In  the  Mediterranean  case  [15],  the  higher  value  pHs8.3  increases  both 
the  Doric  acid  (A2)  and  the  magnesium  carbonate  (A3)  coefficients  Dy  a 
factor  of  two  compared  to  the  N.  Atlantic.  The  curves  do  not  differ  Dy  such 
a  large  factor  at  the  lower  frequencies  because  the  relaxation  frequencies 
are  higher.  The  effect  of  lower  relaxation  frequencies  in  the  Arctic,  for 
the  same  value  pH=8.3,  is  shown  by  the  dsshed  curve. 

The  value  pH-8.0  has  been  assumed  for  Thorp's  experiment  end  is  used  as 
reference  value  in  the  3-relaxaUon  model.  Predictions,  based  on  archival 
pH  values,  show  good  agreement  within  experimental  limits  for  all  regions 
of  the  World  Ocean  examined.  However,  small  adjustments  of  parameters 
can  be  made  if  new  absorption  and/or  pH  data  indicate  the  need. 

Details  of  the  model  In  the  Arctic  regions  will  be  considered  "ext. 
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Arctic  Environment 


Figure  5.  Arctic  Oceon  feotures 


Figure  5  shows  the  general  features  of  the  Arctic  Ocean.  Depths  greater 
than  3  km  occur  in  only  five  main  regions,  namely,  the  Canadian,  Makarov, 
From,  Nansen  and  Norwegian  basins.  In  these  regions,  good  propagation  via 
the  half-channel  can  be  expected  to  very  long  ranges.  Spreading  loss  is 
then  cylindrical,  however,  scattering  losses  limit  useoble  frequencies  to 
less  than  100  Hz  for  ranges  much  greater  than  100  km. 

The  basins  are  separated  by  the  various  ridges,  the  most  prominent  being 
the  Lomonosov,  which  runs  through  the  pole.  Surrounding  the  deep  basins 
are  the  shelf  regions,  in  which  the  depths  rapidly  decrease  to  less  than 
500  m.  In  the  shallower  waters,  the  bottom  can  cause  even  higher  losses 
because  the  deeper-going  rays  ere  cut  off. 
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Arctic  Attenuation 


Frequency  (kHz) 

Figure  6.  Arctic  attenuation  measurements. 

Figure  6  compares  measured  attenuation  dais  with  estimated  absorption 
for  central  Arctic  waters  The  solid  circles  are  1959  data  ( 16j  and  open 
circles  areTristen/Fram-82  data  { 17) 

The  striking  features  of  the  data  8re  the  magnitude  and  resemblence  to 
relax8tional  absorption.  In  the  low-frequency  range,  values  are  more  than 
10  times  greater  than  predicted  absorption  The  relaxation-like  behavior 
is  believed  due  to  experimental  conditions.  Attenuation  usually  limits  the 
useable  range  to  less  than  100  km  at  higher  frequencies  and  absorption 
dominates  in  this  regime.  Smaller  losses  at  low  frequencies  require  much 
longer  ranges  for  measurement  and  scattering  then  becomes  dominant 
The  likely  low-frequency  loss  mechanism  is  scattering  at  the  ice-water 
Interface.  Scattering  theory  shows  that  loss-per-bounce  depends  almost 
linearly  on  grazing  angle  Since  the  skip  distance  has  roughly  the  same 
dependence,  loss-vs-range  tends  to  depend  more  on  gradient  than  angle  for 
grazing  angles  less  than  about  6° 

Larger-angle  rays  encounter  a  much  smaller  gradient,  however,  initial 
divergence- loss  is  greater  and  they  tend  to  contribute  little  at  the  longer 
ranges  Cutoff  by  the  bottom  may  be  an  important  factor 


The  first  attempts  to  model  underice  scattering  were  based  on  the  theory 
of  Marsh  [16].  The  Marsh  formula  for  forward-scatter  loss  employed  the 
Neumann-Pierson  sea-surface  spectrum  and  the  assumption  was  made  that 
the  ice-water  interface  acts  like  a  pressure-release  surface. 

The  theory  appeared  to  give  plausible  results;  i.e.  values  consistent  with 
experiment  were  predicted  for  the  RMS  standard  deviation  2.4  m  for  the 
underice  roughness.  This  value  was  in  fair  agreement  with  measurements 
of  underice  profiles  by  upward-looking  SSN  sonar. 

The  approximations  in  harsh  theory  make  it  valid  over  a  limited  range  of 
frequencies  end  grazing  angles.  Recent  advances  in  scattering  theory  by 
Brekhovskikh  &  Lysenov  [19]  now  permit  calculation  of  forward-scatter 
loss  for  all  angles  and  frequencies.  To  carry  out  the  necessary  numerical 
calculations,  it  is  convenient  to  have  an  analytic  statistical  model  of  the 
underice-roughness. 

In  the  Tristen/Fram-82  experiment  [17],  underice  profiles  were  obtained 
concurrently  with  the  propagation  measurements.  Profile  data  for  400  km 
range  were  analyzed  in  1  km  segments  and  averages  were  taken  over  ten 
contiguous  segments.  The  power  spectra,  correlation  functions  and  PDF's 
(probability-density  function)  of  draft  were  calculated  by  the  FFT  method. 
Standard  deviations  were  derived  from  the  PDF’s  and  correlation  lengths 
were  estimated  by  fitting  with  the  model  correlation  function. 

Average  values  of  correlation  length  and  RMS  standard  deviation  were 
found  to  be  44m  and  2m,  respectively,  and  the  scatter  of  the  10  km  date 
was  reasonably  uniform.  Correlation  between  standard  deviation  end  draft 
was  quite  high;  however  plots  of  correlation  length  vs  standard  deviation 
showed  considerable  scatter. 

Figure  7  is  a  skeich  of  a  small  segment  of  an  underice  profile.  The  power 
spectrum  shows  topical  FFT  data  and  the  curve-fit  with  the  1 -dimensional 
analytic  spectrum.  The  correlation  function  R(r)  and  the  2-dimensional 
spectrum  have  been  derived  analytically  from  the  1 -dimensional  spectrum. 
The  angular-distribution  of  R(r)  has  been  taken  to  be  effectively  isotropic 
and  K, is  the  modified  Bessel  function. 

The  statistical  parameters  are  correlation  length  and  standard  deviation 
The  original  theory  of  Marsh  required  only  standard  deviation  because  the 
correlation  length  is  implicit  in  the  Neumenn-Pierson  spectrum. 


n  9  2  -3/2 

1- Oimensionel  Spectrum  S^h^Q^QQ+Q  ) 

Correlation  Function  R=Q0rKj(Q0r) 

2- Dimensional  Spectrum  S2=  2  h2Ct§  (oj  ♦Q2 )  2 

Standard  Deviation  RMS  h=2m 
Correlation  Length  L=2/Q0=44m 


Figure  7:  Ice  roughness  spectrum  model 


Scattennc-loss  theory 


The  reflection  coefficient  G  is  token  here  as  the  ratio  of  the  coherent 
signal-intensity  in  the  specular  direction  to  the  incident  intensity.  In  the 
perturbation  approximation  we  let  Q=  t -SL  where  SL«1.  By  removing  the 
Neumonn-Pierson  spectrum  from  the  Marsh  formula,  his  result  becomes 
SL=3.3  k3^2  L",/2  h2  sin6,  where  6  is  the  grazing  angle,  ic=2tff/C  is  acoustic 
wavenumber  and  C  is  the  sound-speed 

The  rsy  skip-range  is  given  by  R(B)=2C  tan8/g.  Loss  vs  range  is  therefore 
almost  independent  of  ray-angle  for  constant  gradient  g.  Taking  g=.G6/sec. 
(near  the  surface)  and  the  values  h=2m,  L=44m,  the  attenuation  coefficient 
becomes  A=4.34  SL/R(8)  *1.5  f3/2  dB/km,  where  f  is  in  kHz. 

This  approximation  gives  good  agreement  with  the  data  trend  of  Figure  6 
but  it  is  low  by  roughly  a  factor  of  two.  The  discrepency  is  due  mainly  to 
the  lower  value  of  rms  roughness.  More  exact  computation,  using  computer 
codes  for  estimating  propagation  loss  and  adding  the  calculated  reflection 
loss-per-bounce,  gives  comparable  results. 

The  Marsh  formula  is  valid  only  over  a  limited  domain.  Since  it  involves 
correlation  length,  the  scatter  beam-pattern  is  clearly  involved.  A  general 
theory  for  all  domains  is  given  by  equation  9.6.3  of  B&t,  which  is  a  double 
integral  over  the  surface-roughness  spectrum. 

In  this  paper,  the  B&L  integral  is  transformed  into  scatter-angle  space 
to  sho  w  where  the  energy  goes.  For  the  analytic  ice  spectrum,  the  problem 
then  reduces  to  a  single  integral  to  be  evaluated  numerically.  Evaluation  of 
the  E&L  double-integral  gives  ioentical  results. 

Asymptotic  solutions  of  the  integral  are  8lso  obtained,  which  show  the 
following  limiting  conditions- 

in  the  Eckan  regime  [20j,  the  scattered  energy  tends  toward  specular 
and  the  losses  go  as  f2.  Values  then  become  higher  than  the  Marsh  formula; 
however,  the  method  of  small  perturbations  (MSP)  requires  the  condition 
kh  sm8  «1,  which  is  violated  at  the  higher  frequencies. 

in  the  Marsh  regime,  the  losses  go  os  f3/2  and  the  scattering  beomwidth 
becomes  important.  Some  energy  will  then  escape  and  be  absorbed  in  the 
bottom,  which  could  explain  why  there  is  no  evidence  of  Incoherent  energy 
at  longer  ranges. 

In  the  Rayleigh  regime,  scattering  becomes  diffuse  and  losses  go  os  f4. 
For  the  correlation  length  44m,  low-frequency  roll-off  can  be  expected 
below  about  50  Hz  and  losses  will  then  be  less  than  the  Marsh  formula. 

The  results  are  summarized  in  the  following  equations: 
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Coherent  Energy  Reflectivity  Q=t-SL 


-14?  */2  2  * 

SL=(2ir)  k  h  sin8  {dQ’sin  e‘cos8‘JdHfrdrR(r)Jn(k'r) 

o  -*  o 


k’=k(cos28+Gos2B'-2cos*cosBcos8’)1/2 
k=2iff/C  R(rt=Q^r  KiO^r) 


k2J  rdrR(r)JQ(k,r)=2(1+(k7Q0)2)  2gFt 
0 

First  Integral 

(2ir)  ‘  fdnF,=(  1  +u2)((  1  +u  jV)'s/z  5F. 
-V  1  4 

■ 

UzCk/Q^)2 (cogWosfr)  Vr2(k/Q0J2 COS0COSB’ 

Second  Integral 

“  }  2 
SL=tkh)  si.i0jd9‘sin  8‘cos0'F2 
o 


Final  integral  to  evaluate  numerically 
Eckart  SL=(2kh  sir.B)2  Ru»1  kL»l 

Marsh  SL=3.3(kh)2sirB(kL)1/2  Ru«l  kL»l 
Rayleigh  SLi:2(kh)2(kL)2s;nB/3  kL«l 

U2/fig  Ru=kL  3in2(8/2) 

Asymptotic  Ar^roximations 

The  theory  is  based  on  Bragg-scattering  wherein  k'  defines  a  matching 
between  the  acoustic  (k)  and  the  surface  (Q)  wavenumbers  as  a  function  of 
the  incident  (G)  and  scatter  (8‘)  grazing  angles. 

The  first  integral  over  radius  r  yields  an  analytic  function  'ovolving  the 
ratio  k'/£50.  (JQ  is  the  Bessel  function.) 

The  second  integral  over  azimuthal  angle  b  is  also  analytic,  leaving  only 
an  integral  over  scatter  angle  8'  to  evaluate  numerically.  The  asymptotic 
solutions  provide  a  check.  (Parameter  Ru  refers  to  Rutherford-scattering.) 


Figure  8:  ice  spectrum  Integrands  vs  scatter  angle. 

Integrands  of  the  scattering-loss  equation  are  a  measure  of  scattering 
directivity.  Figure  8  show  plots  of  amplitude  vs  scattering  grazing-angle 
for  an  incident  grazing  angle  8=10°.  The  curves  have  been  normalized  by 
dividing  by  (kh^sinS  in  order  to  restrict  the  range.  At  low  frequencies,  the 
scattering  becomes  diffuse  in  the  helf-spsce,  while,  at  high  frequencies, 
it  tends  to  become  specular. 

The  integrands  can  be  considered  as  the  energy  scattered  from  an  area  of 
radius  equal  to  the  correlation  length.  It  is  therefore  coherent  with  the 
incident  plane-wave.  Summing  over  a  random  ensemble  of  areas  makes  the 
result  incoherent. 

The  integrands  are  similar  to  conventional  beam~p8ttems;  however,  they 
are  a  measure  of  the  flux  in  a  cone  and  ere  not  the  usual  cross-sections. 
The  fraction  of  scattered  energy  remaining  in  the  refractive  regime  can  be 
readily  estimated  from  the  curves.  For  example,  if  15°  rays  hit  the  bottom, 
more  than  1/2  can  be  lost  for  each  surface-  reflection  at  frequencies  less 
than  1  kHz.  After  several  reflections,  the  incoherent  levels  will  converge 
rapidly  and  become  negligible. 
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Figure  9:  Ice  spectrum  scattering  coefficients  vs  frequency. 

In  Figure  9,  values  of  the  integral  obtained  by  numerical  integration  are 
plotted  as  a  function  of  frequency  for  grazing  angles  4-16°,  which  is  the 
range  of  interest.  Note  that  the  curves  have  been  normalized  by  dividing 
the  scatter-loss  by  (kh^sinS.  This  makes  the  Marsh  formula  (dashed  line) 
and  the  Rayleigh  low-frequency  asymptote  independent  of  angle  while  the 
Eckart  high-frequency  asymptotes  become  proportional  to  sin8.  The  Marsh 
formula  becomes  inaccurate  at  the  lower  frequencies  where  the  Rayleigh 
limit  becomes  important  and  also  at  the  higher  frequencies  where  the 
Eckart  limit  becomes  important. 

loss  estimates  below  100  Hz  by  the  Marsh  approximation  have  already 
proven  low  by  more  then  a  factor  of  2.  The  values  obtained  by  integration 
become  even  smaller  at  low  frequencies.  Since  the  conditions  for  the  MSP 
theory  are  apparently  valid  in  this  range,  the  cause  of  disagreement  is  not 
clear.  The  values  do  increase  in  the  Eckart  regime,  but  the  MSP  conditions 
rapidly  become  invalid  here. 

One  possible  additional  loss  mechanism  is  dissipation  within  the  ice; 
however,  very  high  backscatter  levels  ere  also  observed  end  this  suggests 
that  there  is  more  scattering  than  can  be  accounted  for. 
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Attenuation  summaru 


Figure  10:  Attenuation  vs  frequency  models. 

The  attenuation  results  are  summarized  In  Figure  10,  showing  the  two 
scattering  models  and  the  experimental  data. 

The  solid  curve  on  the  left  Is  the  Marsh  formula  in  which  scattering  loss 
increases  as  f3/2  Predictions  are  low  by  more  than  a  factor  of  two  at  the 
lewer  frequencies. 

The  dashed  curve  is  the  numerical-integration  model  covering  ell  three 
regimes.  The  roll-off  below  100  Hz  indicates  a  gradual  transition  to  the 
Rayleigh  regime  In  which  scattering  becomes  diffuse  end  scattering  loss 
approaches  an  f4  asymptote.  The  increase  in  magnitude  above  200  Hz  is  the 
transition  to  the  specular  Eckart-r§gime  where  the  asymptote  is  f2. 

Clearly,  the  more  general  scattering  model  and  more  accurate  knowledge 
of  the  underice  statistics  has  not  improved  agreement  between  experiment 
and  theory.  Without  the  pressure-release  assumption  for  the  ice-water 
interface,  finite  impedance  effects  would  cause  the  magnitudes  to  fall  off 
even  more  rapidly  below  200  Hz.  Air  trapped  near  the  bottom  of  the  ice  is 
a  possible  mechanism  for  effectively  zero  impedance;  however,  further 
experiments  are  required  for  verification  Even  with  this  assumption,  the 
agreement  with  data  is  still  far  from  satisfactory. 
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Figure  1 1:  Backscatter  strength  vs  grazing  angle. 

From  the  theory  of  Marsh  121),  the  ratio  of  the  scattered  to  the  incident 
intensity  is  proportional  to  k4  sin2  8  S^tQ),  where  Q-k'  and  k'  is  the  Bragg 

wavenumber  For  backscattering,  8'=-8  and  the  bock  -'altering  strength  is: 
SS(dB)=  40  log(stn8)  ♦  10  log[k4  S2(Q»  Q~2k  cos8 

Since  the  strength  goes  os  k4  and  the  spectrum  goes  asymptotically  as  Q"4, 
the  result  is  independent  of  frequency  obove  about  50  Hz  and  is  simply: 
SS(dB)~-30  ♦  40  log  (tar.8) 

The  data  points  in  Figure  1 1  are  everege  values  for  the  band  40-500  Hz. 
The  grazing  angles  are  all  close  to  10°  and  no  trend  is  obvious  for  such  a 
limited  range. 

The  Marsh  formula  prediction  falls  well  below  the  data.  The  two-scale 
model  of  Kur’yanov  (221  attempts  to  account  for  the  effects  of  finite  slope 
by  assuming  that  the  scetterers  are  small  in  scale  compared  to  the  scale 
of  t’.ie  gross  roughness.  The  model  of  Greene  and  Stokes  1231  assumes  that 
the  slope  distribution  of  the  gross  roughness  is  non-Gaussian  and  that  the 
backscattering  Is  dominated  by  regions  of  highest  slope,  i.e.  at  pressure 
ridges  Curve  *6&S*  shows  the  prediction  for  the  slope  value  30°. 
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Figure  12:  Backscatter  strength  vs  frequency. 

Figure  12  shows  the  backscettering  data  vs  frequency  for  grazing  angles 
near  10°.  The  'G&S*  prediction  falls  off  below  50  H2  because  the  acoustic 
wavelength  becomes  greater  than  the  horizontal  scale  of  roughness. 

Finite  impedance  at  the  ice-YY8ter  interface  could  also  have  significant 
effect  at  low  frequencies  here;  i.e.  if  the  impedance  discontinuity  begins 
to  disappear  below  200  Hz,  scattering  in  any  direction  will  rapidly  become 
negligible.  Therefore,  without  the  pressure-release  assumption,  matters 
would  be  even  worse. 

It  is  clear  that  surface  Dackscatter-strengths  in  ice-covered  regions  are 
extremely  high.  It  is  also  clear  that  the  theoretical  models  for  both  the 
backscattering  strength  end  forward-scattering  loss  are  not  in  adequate 
agreement  with  experiment,  even  with  the  pressure-release  assumption 
and  composite-roughness  models 

The  only  logical  conclusion  appears  to  be  that  the  scattering  theory  is  et 
fault.  Further  theoretical  investigations  are  obviously  needed  to  see  if  a 
model  can  be  developed  that  accounts  for  the  observed  magnitudes.  Effects 
of  the  ice  Impedance  must  be  Included  8nd  this  could  require  experiments 
to  determine  the  pertinent  physical  properties. 
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Figure  13:  Comparison  of  pH  profiles. 

Figure  13  compares  an  Arctic  pH  profile  with  the  extreme  cases  in  the 
lower  latitudes.  Although  the  pH  values  tend  to  be  fairly  high,  it  is  clear 
that  sound  absorption  has  a  for  less  impor  tant  role  in  the  Arctic  because 
the  excess  attenuation  is  so  great 

In  earlier  reports  (1,21,  the  problem  of  pH  variability  and  prediction  of 
its  effects  on  absorption  was  addressed  From  analysis  of  propagation  data 
from  the  lower  latitudes,  it  was  shown  that  one  method  of  estimating  the 
depth  variations  is  to  use  selected  depth-values  for  various  propagation 
modes;  i.e.  axial  values  for  sound  channels,  surface  values  for  surface 
ducts  and  2  km  values  for  convergence  zones. 

Contours  of  pH  for  the  surface  and  for  depths  0.5  km  and  1  km  from  the 
World  Ocean  Atlas  Vol.  1  (24)  were  used  together  with  contours  for  2  km 
depth  from  Vol.  2  in  developing  appropriate  charts  for  the  three  modes  of 
propagation.  Since  the  Russian  pH  contour  intervals  are  large  to.  1  pH  unit), 
interpolation  was  required  to  achieve  the  desired  accuracy.  Correction  of 
values  to  in-situ  pressure  was  also  required.  The  sound-channel  charts 
were  based  on  the  analysts  of  Russian  pH  contours  by  Lovett  125)  while  the 
CZ  and  surface  charts  were  derived  mainly  from  the  Russian  report. 


In  coses  where  there  is  no  single  clearly  dominant  mode,  effective  loss 
con  be  estimated  by  integroting  over  oil  roy  poths  using  o  pH  profile.  Five 
contour  charts  were  provided  for  the  selected  depths  0,  .5,  1,2  end  4  km 
for  this  purpose.  The  GEOSECS  (26]  data  have  also  been  used  to  derive  the  2 
km  contours  in  ports  of  the  Pacific  Ocean  not  covered  by  Russian  work. 

The  estimation  problem  at  higher  latitudes  differs  in  several  respects 
because  the  thermocline  disappears,  making  refraction  uniformly  upward. 
Water  depths  tend  to  be  smaller  also.  Since  details  of  pH  variability  nearer 
the  surface  are  more  important,  a  change  in  scale  is  indicated.  PH  contours 
in  the  World  Ocean  Atlas:  Arctic  Ocean  Vol.  3  cover  a  depth  range  0-3  km. 
Only  the  range  0-1  km  is  required  for  calculations. 


A=Al(MgS04)+A2(B(OH)3)+A3(MgC03) 

An=(s/35)  anf 2 

fr./(f2+fn> 

B,  =  O.5x.0-D<km)/20 

T/60 

f  j  =  50x10 

f2  =  0.9x1 0T/7° 

.  T/30 

0  2=  0.1XK 

s3=  0.03xK 

f  3  =  4.5x1  0 

The  pH  parameter  K=l(/pH*8)  has  been  substituted  in  the  global  model 
formula  above.  Salinity  dependence  is  teken  as  S/35,  with  the  caveat  that 
errors  may  be  excessive  outside  the  range  30-40  ppt.  Temperature  profiles 
can  be  derived  from  the  SVP  used  in  the  computer  code.  Actually,  they  may 
not  be  required  because  the  range  is  so  small.  Salinity  variations  appear  to 
be  negligible. 

Figure  14  shows  the  rough  outline  of  the  1  km  depth  contour.  Propagation 
can  be  expected  to  be  absorption-limited  for  ranges  up  to  roughly  50  km  in 
this  region.  In  shallower  waters  outside  the  region,  effects  of  the  bottom 
can  become  significant  and  ranges  are  correspondingly  less 

Figures  15  and  16  show  surface  contours  of  the  K-fector  for  winter  and 
summer,  derived  from  pH  contours  in  the  World  Ocean  Atlas  Arctic  Ocean 
The  complexity  of  the  structure  during  warmer  periods  is  probably  due  to 
melting  and  runoff.  This  variability  should  not  effect  absorption  because 
the  deeper  portions  of  the  ray  paths  are  the  most  critical. 

Figures  17-22  show  the  K-factor  contours  for  ell  the  depths  given  in  the 
World  Ocean  Atlas-  Arctic  Ocean  Seasonal  venations  8re  expected  to  be 
much  smaller  at  depths  of  100m  or  greater  and  the  structure  appears  to  be 
relatively  simple 
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16.  Surface  K  contours-summer 
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300m  K  Contours 


Figure  21.2  km  K  contours 
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figure  22  3  km  K  contours 
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Figure  22:  Arctic  Ocean  K-profiles. 

The  K-contour  method  has  proven  to  give  adequate  accuracy  in  modeling 
the  ectual  profiles  throughout  the  World  Ocean.  The  algorithm  to  be  used 
again  here  was  suggested  by  Dr.  A  H  Nuttall. 

The  K  values  for  D=0,  0.1,  0.3,  0  5,1  km  are  labeled  Dn  where  n=0,l, 2,3,4 

The  profile  is  generated  from  the  equation 

K(D)=K(D4)*  IC0  ♦  C,  D  ♦  C2  D2  ♦  C,  D3  ♦  C4  D4]  expKaD)5] 

where  e=4/fcm  and  D=1.5  are  found  to  give  'best'  results  for  Arctic  regions. 
The  five  equations  to  be  solved  for  the  coefficients  Cn  are  then  given  by: 

V  C1  Dn  *  C2  2  *  ^  Dn  3  *  C4  4  =  "  «D4)1  eXP((4Dn)1  5} 

Either  algebraic  or  matrix  methods  can  be  used. 

Figure  22  shows  typical  profiles  derived  from  the  K-contour  charts  for 
several  regions  in  the  Arctic  Ocean.  The  profiles  appear  quite  consistent 
with  the  Russian  data.  Generally,  the  pH  value  is  greatest  at  the  surface 
and  becomes  constant  a'  depths  greater  than  about  0.5  km  with  a  minimum 
near  0.1  km  depth  in  some  cases.  The  range  of  variation  tends  to  be  small 
compared  to  the  World  Ocean  at  lower  latitudes  1 1 ,21. 
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Figure  23:  Absorption  spectrum  range  in  the  Arctic. 

Figure  23  shows  the  predicted  range  of  variation  of  absorption  spectra 
In  the  Arctic  for  the  pH  range  6.0-6.3  or  K-fector  range  1  -2.  Temperature 
Is  taken  as  -  PC  and  salinity  as  35  ppt.  Effects  of  pH  below  10  kH2  are  due 
mainly  to  changes  In  the  boric  acid  coefficient.  At  very  low  frequencies, 
the  max/mtn  ratio  approaches  the  factor  two.  When  compared  to  Atlantic 
spectra  for  pH=B.O,  the  range  can  exceed  two  because  the  temperature  is 
higher 

The  d8t8  points  are  from  the  Baffin  Bay  experiment  of  OREA  (12  Vol  ill. 
These  are  the  only  available  low-frequency  Arctic  data  that  clearly  show 
the  effects  of  the  absorption  limit  The  experiment  was  carried  out  in  the 
late  summer  under  ice-free  conditions.  Near-surface  temperatures  were 
sufficiently  high  to  form  a  weak  thermocline  and  the  propagation  mode 
was  sound-channel,  the  axis  being  near  100m  in  depth.  The  experimental 
range  extended  to  roughly  400  km.  The  dashed  curve  for  pH=B.2  (Ks  l  .6) 
gives  the  ‘best'  data-fit  at  the  higher  frequencies 

The  excess  loss  at  the  lower  frequencies  is  probably  due  to  scattering  by 
temperature  inhomogenettles  in  the  medium,  since  neither  ice  or  surface 
waves  can  be  held  responsible 
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Model  prediction  comparison 
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Figure  24:  Model  comparison  of  2.5-4  5  kHz  two-way  losses  at  40  km. 

in  most  cases,  propagation  In  the  high  Arctic  will  be  absorption-limited 
only  for  single  refraction  paths.  The  ranges  of  concern  are  therefore  much 
shorter.  Two-way  loss  predictions  for  the  K  model  and  the  Thorp  equation 
for  the  range  40  km  are  compared  in  Figure  24.  The  differences  between 
the  models,  A=A(Kmod)-A(Thorp),  also  include  the  expected  +  dB  error  of 
the  K  model,  which  indicates  the  degree  of  significance. 

The  ray-integretion  method  8na  the  algorithm  profile  were  used  in  the  K 
model  calculations.  Values  of  K  at  the  five  depths  were  fitted  to  generate 
the  profiles  Salinity  was  taken  as  35  ppt  and  temperature  as  -rc.  Losses 
at  selected  frequencies  F (kHz)  were  calculated  using  the  ray- integral  1  on 
method  Relative  errors  were  calculated  concurrently  using  the  expected- 
value  AK=i0.05.  The  Thorp  values  were  also  calculated  concurrently  so  as 
to  minimize  relative  errors. 

Figure  25  shows  similar  calculations  for  frequencies  5-9  kHz  at  10  km 
range  and  for  frequencies  10-15  kHz  at  5  km  range  Differences  between 
models  become  negligible  at  the  higher  frequencies.  Results  above  10  kHz 
are  in  very  good  agreement  with  reported  experimental  values  1271 
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Figure  25  Model  comparison  of  5- 16  kHz  two-w8y  ’osses  at  5  end  10  km 
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Figure  26:  Model  comparison  of  0.5-2  kHz  one-way  losses  at  200  km. 

Ice-free  conditions  can  exist  in  the  Norwegian  Sea  during  the  summer 
months  and  propagation  to  very  long  ranges  is  possible  at  low  sea-states 
if  the  surface  temperature  rises  enough  to  form  a  shallow  sound-channel, 
refraction  effects  can  make  surface  scattering  negligible.  However,  the 
thermocline  structure  can  still  be  ‘rough*  and  this  evidently  produces 
internal  scattering.  Excess  losses  may  tnen  occur  a  lower  frequencies,  like 
that  observed  in  the  Baffin  Bay  experiment 

Figure  26  compares  predicted  one-way  absorption  loss  for  sound-channel 
propagation  in  the  Norwegian  Sea  at  200  km  range  with  other  regions  of 
the  World  Ocean.  The  full  range  of  ph  end  temperature  effects  is  realized 
only  in  these  lower  frequencies. 

From  the  earlier  analysis  (1,21,  absolute  error  of  the  absorption  model  ic 
estimated  to  be  less  than  ±15£,  i.e  RMS  error  of  the  coefficient  in  dB/km 
is  not  expected  to  be  greater  than  this  if  not  limited  Pu  the  accuracy  of 
the  environmental  factors.  Relative  errors  in  estimating  the  effects  of  pH 
and  temperature  can  be  expected  to  be  much  smoiler 
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Conclusion 


The  purpose  of  this  effort  has  been  to  extend  the  global  model  for  sound 
•absorption  in  sea  water  [1,2]  to  cover  Arctic  regions.  The  method  used  in 
estimating  sound  absorption  at  lower  latitudes  appears  equally  effective 
in  polar  waters.  The  only  modification  required  is  the  change  in  vertical 
scale.  Since  thermocline  effects  are  absent,  upward  refraction  dominates 
propagation  and  variations  nearer  the  surface  become  more  important.  The 
depth  range  0- 1  km  has  therefore  been  selected  for  the  K-contours  and  the 
scale  of  the  K-prof ite  algorithm  changed  from  1/km  to  4/km. 

The  recommended  method  is  numerical  integration  of  losses  over  all  ray 
paths  by  computer  methods  employing  K-profiles.  The  five  points  from  the 
K-contour  charts  should  provide  the  required  accuracy  for  the  K-profiles. 
The  algorithm  method  of  generation  is  also  recommended. 

Temperature  profiles  can  likewise  be  generated  from  the  SVP  used  in  the 
computer  code.  Since  the  nominal  range  is  only  -2&C  to  +4°C,  effects  will 
be  small.  The  default  value  -  1°C  is  suggested 
Salinity  variations  fall  in  the  range  32-35  ppt  (except  near  the  surface) 
and  effects  will  also  be  minimal.  The  default  value  35  ppt  is  suggested  if 
local  data  are  unavailable 

The  K-profiles  can,  of  course,  be  generated  graphically,  however,  making 
a  fit  involves  subjectivity  and  translation  to  a  computer  code  is  difficult. 
Linear  interpolation  would  probably  be  more  effective 
A  much  simpler  but  less  accurate  method  of  estimating  absorption  is  to 
select  one  K-vnlue  from  the  depth  contour  best  suited  to  the  propagation 
conditions  in  question.  In  deeper  water,  for  example,  the  0  5  km  contours 
would  probably  be  sufficiently  accurate  While  the  method  is  useful  for 
rapid  calculation,  it  is  also  more  subjective 
The  final  caveat  is  that  no  reliable  date  for  polar  waters  appear  to  be 
available  except  for  frequencies  above  10  kHz  Absorption  predictions  at 
the  lower  frequencies  therefore  rely  solely  on  extrapolation  of  results 
obtained  at  higher  temperatures  Single-refraction  path  experiments  in  the 
range  1-10  kHz  would  help  to  resolve  this  problem.  Excess  losses  observed 
at  tower  frequencies  could  involve  some  unknown  absorption  mechanism 
that  becomes  evident  only  at  near-freezing  temperatures  Experiments 
using  the  resonator  technique  would  serve  os  a  method  of  investigating 
this  possibility 
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